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We study the static, spherically symmetric black hole solutions for a non-minimally coupled multi- 
scalar theory. We find numerical solutions for values of the scalar fields when a certain constraint 
on the maximal charge is satisfied. Beyond this constraint no black hole solutions exist. This 
constraint therefore corresponds to extremal solutions, however, this does not match the n = 
constraint which typically indicates extremal solutions in other models. This implies that the set 
of extremal solutions have non-zero, finite and varying surface gravity. These solutions also violate 
the no-hair theorems for N > 1 scalar fields and have previously been proven to be linearly stable. 
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I. INTRODUCTION 

Black hole solutions with scalar fields are usually con- 
strained to possessing only secondary hair by the "no- 
hair conjectures" Early attempts to find black hole 
solutions coupled solely to a scalar field found solutions 
where the scalar field diverged at the putative horizon Q . 
Hence, technically such solutions cannot actually be con- 
sidered black hole solutions due to a lack of a regular 
horizon. They are unlikely to have any physical rele- 
vance [3j. The original "no- hair conjectures" have since 
been violated in a number of cases, either via coupling the 
scalar fields to both gravity and gauge fields, or through 
violation of the dominant energy condition. When the ex- 
istence of scalar hair depends on a non-vanishing gauge 
field, and is entirely fixed by the mass, gauge charge and 
angular momentum this is called secondary hair [J]. In 
this paper we discuss a solution with contingent primary 
hair [f| , that is to say the scalar hair depends on the ex- 
istence of a non- vanishing gauge field but its behaviour is 
not entirely fixed by the values of the other asymptotic 
parameters. We briefly list some further "hairy" black 
hole solutions. 

• Scalar fields coupled to higher order gravity have 
been heavily investigated since they arise naturally 
in low energy effective 4-dimensional string theo- 
ries. The Gauss-Bonnet term, which is the only 
ghost-free leading order curvature correction, has 
naturally been of particular interest @, Q ■ 
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• Minimally coupled scalar fields with dominant en- 
ergy condition violating potentials have been shown 
to allow non-trivial hair 

001, EH, E3- Examples 
have been found both analytically and numerically 
provided there is at least one global minimum with 
V(4>) < 0. 

• Theories which couple gravity to non-Abclian 
gauge fields such as Einstein- Yang-Mills, Einstein- 
Yang- Mills-Higgs and Einstein-Skyrme, usually 
contain nonlinear self-interactions and admit 
"hairy" black holes. Einstein- Yang-Mills-Higgs and 
Einstein-Skyrme also include scalar fields. These 
vanish exponentially at infinity, however, and thus 
they do not have "Gauss-like" scalar charge. These 
hairy black holes were thought to be generally un- 
stable but it has been shown that some branches 
of solutions of the Einstein-Skyrme black holes are 
linearly stable [6]. Whether they are non- linearly 
stable remains an open question. 

• Scalar fields non-minimally coupled to an Abclian 
gauge theory have been shown to have hairy solu- 
tions [H, 0, [ID, HH . Such theories arise naturally 
in Kaluza-Klein theories and effective low-energy 
limits of string theory with a non-trivial dilaton. 

Despite the solutions listed above being beyond the 
premises of the original "no-hair conjecture" , they are 
still considered interesting as tests of the limits of the 
conjectures. Stability is still an open problem in most 
cases. 

For a non-rotating, static black hole with a single scalar 
field coupled to the U(l) electromagnetic gauge field [HI]. 
This solution is not a member of the Reissner-Nordstrom 
class but is entirely specified by the values of M, Q and P. 
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Adding, an extra scalar field was shown to give more free- 
dom [17J and a version of scalar hair that falls between 
the definitions of primary and secondary hair. This was 
called contingent primary hair and has been generalised 
to N scalar fields with linear stability being shown Q. 
Here we present numerical solutions to this model and 
discuss some of the features. As we will see, the most in- 
teresting feature is that these solutions limit to non-zero, 
finite surface gravity for an extremal black hole solution 
with a general coupling. 

This paper is organised as follows. Section HT1 discusses 
the model to be used and the analytical constraints that 
can be placed on the solutions. Section ITTT1 contains the 
main results and gives details of the solutions found while 
the thermodynamic behaviour of these solutions is shown 
in section IIVI We conclude with a discussion in section 
IVl We use the notation of 5] and define c = Attk 2 = 1. 



solutions to exist we take the magnetic monopole field 
ansatz 



Psin6> 



(4) 



where P is the magnetic charge. This choice is made out 
of convenience. Due to the scalar coupling to the elec- 
tromagnetic sector, the electric ansatz includes depen- 
dence on the scalar fields and is therefore non-trivial, the 
magnetic ansatz, being the 9(f> components of the electro- 
magnetic tensor, avoids these complications. There is no 
longer a simple duality between the magnetic and elec- 
tric solutions although solutions for the electric solution 
should still be tractable if the magnetic solutions exist. 
We could, of course, also consider a situation where both 
are non-zero. As in the single scalar field of [l5[ the scalar 
fields will necessarily vanish if Q = P = 0. 

The G t ' component of the Einstein equations gives 



II. MODEL 
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The general Lagrangian density for the iV-scalar field 
case is 
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where R is the Ricci scalar and is the electromagnetic 
field strength. Initially we consider the N = 2 case with 
no cosmological constant, ie. A = 0. For simplicity we 
split the representation of the scalar fields such that 
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Since there is no potential dependent on any of the scalar 
fields, the Lagrangian density has the same scale invari- 
ance as the Gibbons-Maeda solution [l5| . This invariance 
applies under global re-scalings of the metric g a b w 2 g a b 
where V a u> = 0. 

We use a standard metric ansatz for static, spherically 
symmetric jSchwarzschild coordinates following the for- 
malism of . 

ds 2 = -A 2 (r) fl-^V 



Aje-agi* + X 2 e~ 2 ^ \ P 2 

X 2 + X 2 )r*- (5) 



The linear combination G t * — GJ' of components of the 
Einstein equations gives 

A 



A 
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while the two scalar field equations are 



d r 
and 



1 - 
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Ar'dr* = —AXfgie^ 91 ~y (7) 



d r ((l- ^-^j Ar 2 d r *) = -AX 2 g 2 e- 2 ^^. (8) 

We note these generalise to N + 2 field equations for the 
system given in ([I}, 
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= 1 
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where dfl 2 - Aa2 — 2 
Misner-Sharp mass function. In order for non-trivial 



d6* 2 + sin^ 9d(j} 2 and m(r) is the familiar 
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FIG. 1: Solution given by Ai = A2 = gi — gi = 1.0, P — 0.5 and with horizon values m(rh) = 0.25, $(r/j) = 1.0 and 
9(r h ) = 0.1. Pcx) = 0.8611 and Madm = 0.6674 while A(r h ) = 0.6490 as found by the shooting method: (a) m(r); (b) A(r); 
(c) $(r); (d) *(r). 



III. NUMERICAL SOLUTIONS 

The N = 2 solutions can be obtained numerically with 
the help of the following expansions near the horizon; 

m(r) = m/j + mi(r - r;j) + m 2 (r - r/j) 2 + 

A(r) = A h +A 1 {r-r h ) + A 2 {r-r h ) 2 + ..., 

$(r) = $h. + $i(r -r h ) + <S> 2 {r - r h f + 

*(r) = ^ h + m 1 {r-r h ) + ^ 2 {r-r h f + .... (12) 

In Fig. [1] we show the general form of the solutions for a 
non-extremal case. 

The solution is uniquely defined by 3 asymptotic 
charges in the N — 2 case and by N + 1 charges in the 
general case. The Arnowitt, Deser and Misner (ADM) 
mass, Madm , is given by the asymptotic value of m(r), 
while the asymptotic Gauss-like magnetic charge is given 
by 



Poo = 



(Af + A!) 



(13) 



These two asymptotic charges along with the coefficient 
<j>_i of the 1/r term in the asymptotic $ expansion 



<f>_l <f>_2 



uniquely define the solution. As shown in is 
constrained in the N = 2 case by 



+ tfij + 2M 
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and in general by 

AT JV 



i=l 



i=l 
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(15) 

where 3>_i,j denotes the 1/r coefficient of the ith scalar 
field. This constraint limits the system to N + 1 degrees 
of freedom, with these being in the magnetic monopole 
case, M ADM , Poo and {$-1,1, $_ 1)2 ,...,*-i,jv-i}. The 
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constraint, (|15|) . holds throughout our numerical work to 
the numerical accuracies required. 

We use the shooting method to find solutions such that 



lim A = l, 



(16) 



where this condition is adhered to with a numerical ac- 
curacy of 10 -7 . This exploits the rescaling freedom in 
A(rh) which can be seen in (O-©. This is a necessary 
requirement for the time coordinate t in ([3]) to correspond 
to the proper time of static observers at infinity and gives 
the correct normalisation of the surface gravity. The nu- 
merical limits used to define the asymptotic region are 
m'(r) < 10~ 8 and r > 500r h . 

We have also found solutions to both the iV = 3, A = 
and N = 2, A = — 1 cases. The N = 3 case appear 
elsewhere [19| as there are no particular additional fea- 
tures when compared to the solutions in fig. [TJ From 
these results, however, we would assume that solutions, 
with N > 2 scalar fields, exist, having N + 1 degrees of 
freedom. This may be of interest to string theory mo- 
tivated work, where, in many cases a large or infinite 
number of scalar fields appear in the low energy effective 
4-dimensional theory (see [2d ] for a review). 

The anti-de-Sitter (adS) solutions, while attainable, 
have distinct numerical issues related with finding solu- 
tions over a wide parameter range. However, the critical 
temperature generally exhibited by adS solutions, due 
to the thermal bath, may result in interesting behaviour 
when Hawking evaporation is considered if the unique 
thermodynamic features of the N — 2, A = system 
shown below are also manifest in the A = 1 case. This is 
left to further work although again solutions for a given 
parameter range appear in [l9j . 



IV. THERMODYNAMIC BEHAVIOUR 



The surface gravit 
dinates is given by [18 1 



22j for a black hole in these coor- 
I 

A{r h ) 



4m(rh) 
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Clearly we will have zero-temperature black hole solu- 
tions (k = 0) if m'(r/j) = 1/2. In the horizon expansion 
given above this would correspond to mi = 1/2 while 
from the equation of motion ([5]) we find 



P 1 



mi 



2r 2 

Hence, k = when 



m(r h ) 



s=0 



2\~ 



(18) 



(19) 



where we have used r h — 2m(r h ). We have denoted this 
limiting case as m(r/ l )| ts= o as it turns out not to be the 



extremal case. The "extremal solution" is the solution 
existing with the maximum electromagnetic charge for a 
given mass and scalar charge. 

We find that the separate condition given in Mignemi 
and Wiltshire @], 



N 



PL < M 



ADM 



E* 2 



(20) 



is the constraint for extremal black holes. This con- 
straint has no thermodynamic significance but the equal- 
ity indicates the degenerate horizon. This is a novel 
situation as when (TH))) and (|2"0)) are considered we find 
Poo, extremal < Poa,K=o- For the degenerate horizon, given 
by the equality in ([2"0")l . we find that the horizon becomes 
singular and hence does not have a well-defined surface 
gravity. This is indicated by divergence of the coefficient 
mi in (TT2"]) , 

The limiting cases, however, still indicate extremal 
black hole solutions with non-zero, finite surface gravity 
for general couplings, gi. Such results have been found 
previously for specific couplings [TEl ]. 

The limiting behaviour due to (|2"0|) is shown in fig. [21 
where <I>_i is defined by 



-Adr, 



(21) 



and similarly for 

Contour plots of the surface gravity show Reissner- 
Nordstrom-like solutions when gi = g 2 = 1.0 and Kaluza- 
Klein-like solutions when g\ — g 2 — 3.0 see Fig. [2Ja) 
and Fig. [3jb) respectively. The 'specific heat', de- 
fined as C = (dM ADM /dT)\ Poo ~ (dM ADM /dn)\ P ^, 
changes sign in fig. [21a) while it is always negative for 
fig. (21b). Unfortunately we do not possess the compu- 
tational power to find the limit of the coupling gradients 
that produce these two types of solutions. 

The extremal limit in these cases is not an 'isotherm' 
but instead tends to finite non-zero values where the sur- 
face gravity is decreasing with increasing P^. The con- 
tours mimic those found in [l5j but with a region ex- 
cluded due to the constraint ([2H)l . Fig. 2] shows this 
graphically. However, we would caution that the ex- 
tremal solution falls into a different class from those 
solved by the numerical method implemented here. As 
the value of m'(r^) diverges we do not have well-defined 
solutions and hence no surface gravity. We therefore 
only comment on the limiting behaviour for the extremal 
cases. 



V. DISCUSSION 

We have numerically demonstrated linearly stable 
black hole solutions with contingent primary hair. The 
condition (Tr5|) , as previously derived in [f| , gives N + 1 
asymptotic charges for N scalar fields with two being 
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FIG. 2: Behaviour of constraint given in (|20[) for Ai = A2 = gi — g2 — 1.0: (a) Poo = 1.2; shows the solution approaching 
the limit as it reaches the extremal case at Madm = 0.849 (numerical accuracy of 10~ 5 used for finding Poo and Madm), 
note that although (a) shows the expected behaviour it does not allow comparison with the expected values for k = as all 
quantities are defined at infinity; (b) P = 0.5; approaches the limit m(rh) = 0.233 while the required value for k — from (|19[) 
is m(r h ) = 0.173. 




FIG. 3: Contour plots showing the behaviour of the surface gravity; Ai = X2 — 1.0 with horizon values 4>(rft) = 1.0 and 
^(rh) = 0.1: (a) g± = g^ = 1.0; (b) g± = g^ = 3.0. Note the change in behaviour from Reissner-Nordstrom-like solutions in (a) 
where the specific heat changes sign for Poo = const as it moves away from the extremal limit while in (b) the specific heat for 
Poo = const is always negative and mimics the well-known Kaluza-Klein examples. The extremal limit is shown as a dotted 
line as it is not an 'isotherm', instead the surface gravity decreases with 



Madm and Poo (in the non-zero magnetic monopole so- 
lution considered here) with the other iV— 1 charges being 
the 1/r coefficients of the asymptotic expansion of N — 1 
of the scalar fields, {$_i i i,$_i i 2,...,$_i,jv-i}- This vi- 
olation of the no-hair theorems is, however, not entirely 
within the confines of the premise under which the theo- 
rems were originally derived as we have non-minimal cou- 
pling between the scalar field and the £7(1) gauge field. 

The solutions here may help to shed some light 
on black hole solutions to the low energy effective 4- 
dimensional string theory when coupled with further cor- 
rections, such as, higher-order gravity terms [H, another 
£7(1) field or the inclusion of scalar potentials. Chen et 
al. have found constraints on the value of the coupling in 



the single scalar case when a Gauss-Bonnet term is intro- 
duced. This appears to limit the applicability in string 
theory motivated situations. However, this constraint 
would possibly be weakened by additional scalar fields, 
similar to the case for the slope of the potential when 
additional scalar fields are considered in cosmology [2l| . 

The result of major interest is that the solutions are 
bounded by (j2"0|) and do not contain the K = case. 
At the extremal limit no surface gravity is defined, it 
does, however, limit to finite, non-zero values. Previously 
this behaviour has only been seen in [l5| when consider- 
ing the limiting case between Reissner-Nordstrom and 
Kaluza-Klein type solutions for a single scalar field with 
g = \J D — 1. Although (fTS)) limits the number of inde- 
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pendent asymptotic charges, it does not allow us further 




insight into the nature of the horizon. The Gibbons- 
Maeda solution allowed analysis of the horizon in the 
extremal case, indicating a singularity. This would seem 
likely in the present case as m'(r) diverges at the horizon 
indicating a curvature singularity. 
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FIG. 4: The dotted line shows the constraint (|20|l while the 
solid line indicates the limit imposed by (|19|l for Ai = A2 = 
g\ — g2 — 1.0, $(r"h) = 1.0 and 9(rh) = 0.1 while varying 
over horizon values for m(rh) and the 'bare' P values. 
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